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I. INTRODUCTION 



Measurements of parity nonconservation (PNC) effects in intermediate and lieavy atomic 
systems might provide very stringent tests of the Standard model of electroweak interactions. 
Nuclei and electrons are bound into atoms by electromagnetic interactions that do not 
violate parity. Parity nonconservation arises from the interference of the 7 and gauge 
boson exchange between the nucleus and the atomic electrons. The dominant part of the 
PNC interaction results from the coupling of the axial electronic current to the vector nuclear 
current, and it allows normally forbidden transitions. The nuclear vector current is conserved 
and the PNC effects provide a measure of the electroweak coupling of the elementary quarks. 
An extensive analysis of the sensitivity of a variety of low-energy PNC observables to possible 
new physics phenomena has been recently reported in Ref. [|l|. 

The nucleus acts as a source of the weak charge Qw The quantity which is measured in 
experiments is the electric dipole amplitude between two electronic states which, in the ab- 
sence of the PNC electron-nucleus interaction, would have the same parity. This observable 
can be parameterized [^0] 

GpNc = iQw = i Qw^°'^ + ^Qw ^ ' (1) 

where ^ is a coefficient that depends on atomic structure, Q^'^^"'^ is the Standard model 
prediction for the weak charge, and AQ^^ results from the difference between the neutron 
and proton density distributions in the nucleus. At tree level 

Q^^-^'°''={l-Asm^dw)Z-N. (2) 

Z and N are the numbers of protons and neutrons, respectively, and 6y/ is the weak mixing 
angle. The finite size of the nucleus modifies Z — > qpZ and — > g„A^, where 

qn{p) = J f{r)pn{p){r)£r. (3) 

The function /(r) describes the spatial variation of the electron axial transition matrix 
element over the nucleus. The difference between the neutron p„ and proton pp density 
distributions modifies the effective weak charge 



^Q-^^ = N{l-qJq,). (4) 

The two principal factors, therefore, which hmit the precision with which the weak charge 
Qw can be determined from atomic experiments are: (i) the accuracy of atomic theory cal- 
culations in Eq. ([l|) and the function /(r) in Eq. (|^)), and (ii) the nuclear structure effects, 
i.e. the uncertainties in the neutron and proton density distributions. At present, the level 
of accuracy of atomic structure calculations constitutes the predominant uncertainty in the 
interpretation of atomic PNC observables [0]. For example, the most accurate experimental 
result for the weak charge of the Cs nucleus is Q^^* = —72.QQ{2S)expt{'i^)theory 0- The 
Standard model value including radiative corrections is Qy^'^'^°'^ = —73.20(13). It has been 
suggested that the measurement of PNC effects along an isotope chain might represent a 
useful method to circumvent the atomic theory uncertainty. The ratios of the electric dipole 
amplitudes for different isotopes should not depend on the details of atomic structure. In 
addition to the obvious N dependence, these ratios will be sensitive to the value of the 
Weinberg angle and, eventually, to additional parameters of a more general electroweak the- 
ory [|l],^ . The shortcoming of this method, however, is that the ratios of PNC observables 
will also be sensitive to changes in the neutron distribution p„ along the isotope chain. The 
PNC interaction is directly proportional to the overlap between electrons and neutrons. And 
while extremely accurate data on proton distributions in nuclei are obtained from electron 
scattering experiments, data of comparable precision on neutron density distributions are 
presently not available. It is much more difficult to measure the distribution of neutrons, 
though more recently accurate data on differences in radii of the neutron and proton density 
distributions have been obtained ||§J^. 

A number of studies have been reported on the nuclear structure effects in the inter- 
pretation of atomic PNC experiments in Cs and Pb isotope chains. The basic theoretical 
framework has been analyzed in Ref. 0, and it was pointed out that the uncertainties in 
the neutron distributions would have to be significantly reduced, if isotope ratios are to be 
used in high-precision evaluation of the atomic PNC parameters. On the other hand, it was 
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also emphasized that atomic PNC experiments might provide a model independent tool for 
studying variations in neutron distributions. In a more extensive analysis of electroweak pa- 
rameters and nuclear structure effects , a number of relativistic and nonrelativistic nuclear 
model calculations were performed for the chain of Pb isotopes. The conclusion was that the 
spread among model predictions for the sizes of nuclear structure effects may preclude using 
Pb isotope ratios to test the Standard model at better than a 1% level. The theoretical anal- 
ysis included the Hartree-Fock-Bogoliubov model with Gogny forces, various Hartree-Fock 
Skyrme effective interactions, as well as the relativistic mean-field model. The Hartree-Fock 
Skyrme model, with the inclusion of deformation effects, was used in Ref. to calculate the 
proton and neutron densities in the isotope chain ^^^Cs - ^^^Cs. The good agreement with the 
experimental charge radii, binding energies, and ground-state spins, led to the estimate that 
the uncertainties in the differences of the neutron radii may cause uncertainties in the ratio 
of weak charges of less than 10~^, safely smaller than the anticipated experimental error. 
In a detailed analysis of Ref. [|10| it was shown that extremely high-precision measurements 
of PNC effects in atomic Yb appear possible, and it was suggested that a comparison of 
PNC observables in the wide range of stable isotopes of Yb may provide a unique test of the 
Standard model. The recent experimental result for the weak charge of Cs has prompted 
the analysis of the effects of neutron spatial distribution on atomic PNC in Cs [H. Several 
theoretical models have been used to calculate the modifications to the nuclear weak charge 
due to small differences between the spatial distributions of neutrons and protons in the Cs 
nucleus. It has been shown that the nuclear structure uncertainties are smaller than those 
arising from atomic theory calculations, but are comparable to the present uncertainties in 
the Standard model prediction. The sensitivity of atomic PNC to nuclear structure effects 



in single isotopes and isotope chains has also been extensively discussed in Refs. |1| and [jTT 

In the comparison of the theoretical models that have been used in the analysis of nuclear 
structure effects, it has been pointed out that relativistic mean field models typically 
generate significantly larger neutron radii, and thus would predict much larger corrections 
to the nuclear weak charge. Models of nuclear structure based on relativistic hadrodynam- 



ics [|T2[ have attracted much interest in the last decade. In particular, the relativistic mean- 
field model has been very successfully applied in the description of a large variety of nuclear 
structure phenomena in spherical and deformed nuclei |T3|. The model has more recently 



been extended to include both mean-field and pairing correlations in a unified framework: 
the relativistic Hartree-Bogoliubov model [ITlllTS. A unified and self-consistent description 



of the long range mean-field correlations and the pairing field is especially important for the 
structure of open shell nuclei further away from the line of /5-stability. In addition, in the 
last few years several new and more accurate parameter sets of meson masses and meson- 
nucleon coupling constants for the effective mean-field Lagrangian have been derived. These 
new effective forces have been used not only for a more accurate description of properties 
of /3-stable nuclei, but they have also been applied in studies of isovector properties of ex- 
otic nuclei with large isospin values, of the structure of superheavy nuclei, and properties 
of excited states in the relativistic random-phase approximation and the time-dependent 
relativistic mean-field model. 

In a recent study |TB[ we have applied the relativistic Hartree-Bogoliubov model to the 



description of parity- violating elastic electron scattering (PVES) on neutron-rich Ne, Na, Ni 
and Sn isotopes. For the elastic scattering of 850 MeV electrons on these nuclei, the parity- 
violating asymmetry parameters have been calculated using a relativistic optical model with 
inclusion of Coulomb distortion effects. The asymmetry parameters for chains of isotopes 
have been compared, and their relation to the Fourier transforms of neutron densities has 
been studied. The important impact that PVES measurement of neutron densities could 



have on atomic PNC experiments has been recently discussed in Ref . [|TT| . Since the sensitiv- 
ity of the neutron distribution appears to be approximately the same for PVES and atomic 
PNC, it would of course be ideal if density distributions for PNC experiments were deter- 
mined from PVES. In the present work we use the relativistic Hartree-Bogoliubov model 
to calculate ground-state proton and neutron distributions for the isotopic chains that have 
been suggested for atomic PNC experiments: Cs, Ba, Yb and Pb. Calculated quantities 
will be compared with available experimental data, as well as with previous analyses of 



the nuclear structure corrections to the weak charge In Section II we present an 

outline of the relativistic Hartree-Bogoliubov model. The neutron densities for atomic PNC 
experiments are calculated and analyzed in Section III. 



II. THE RELATIVISTIC HARTREE-BOGOLIUBOV MODEL 

The relativistic mean field theory is based on simple concepts: nucleons are described 
as point particles, the theory is fully Lorentz invariant, the nucleons move independently in 
mean fields which originate from the nucleon-nucleon interaction. Conditions of causality 
and Lorentz invariance impose that the interaction is mediated by the exchange of point-like 
effective mesons, which couple to the nucleons at local vertices. The single-nucleon dynamics 
is described by the Dirac equation 

a ■ V + /3(m + g,a) + g^uo"" + gpT-^p^ + e ^^ 1 = e,i)i. (5) 
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0", oj, and p are the meson fields, and A denotes the electromagnetic potential, g^ g^j, and gp 
are the corresponding coupling constants for the mesons to the nucleon. The lowest order 
of the quantum field theory is the mean-field approximation: the meson field operators are 
replaced by their expectation values. The sources of the meson fields are defined by the 
nucleon densities and currents. The ground state of a nucleus is described by the stationary 
self-consistent solution of the coupled system of the Dirac (|^)and Klein-Gordon equations: 

-A + ml] a(r) = -g^ p,(r) - g2 a\r) - g^ a'(r) (6) 

-A + ml]uj\r)=g^p,{r) (7) 

-A + m^] p%r)=gpp,ir) (8) 

-AA°(r) = ep,(r), (9) 

for the sigma meson, omega meson, rho meson and photon field, respectively. Due to charge 
conservation, only the 3rd-component of the isovector rho meson contributes. The source 
terms in equations (^) to are sums of bilinear products of baryon amplitudes, and they 
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are calculated in the no-sea approximation, i.e. the Dirac sea of negative energy states does 
not contribute to the nucleon densities and currents. Due to time reversal invariance, there 
are no currents in the static solution for an even-even system, and therefore the spatial 
vector components oj, ps and A of the vector meson fields vanish. The quartic potential 

U{a) = ^nila^ + ^g2a' + ^g3a^ (10) 

introduces an effective density dependence. The non-linear self-interaction of the a field is 
essential for a quantitative description of properties of finite nuclei. 

In addition to the self-consistent mean-field potential, pairing correlations have to be 
included in order to describe ground-state properties of open-shell nuclei. For nuclei close to 
the /3-stability line, pairing has been included in the relativistic mean-field model in the form 
of a simple BCS approximation For more exotic nuclei further away from the stability 



line, however, the BCS model presents only a poor approximation. In particular, in order 
to correctly reproduce density distributions in neutron-rich nuclei, mean-field and pairing 
correlations have to be described in a unified framework: the Hartree-Fock-Bogoliubov model 
or the relativistic Hartree-Bogoliubov (RHB) model. In the unified framework the ground 
state of a nucleus |$ > is represented by the product of independent single-quasiparticle 
states. These states are eigenvectors of the generalized single-nucleon Hamiltonian which 
contains two average potentials: the self-consistent mean-field T which encloses all the long 
range particle-hole (ph) correlations, and a pairing field A which sums up the particle- 
particle (pp) correlations. In the Hartree approximation for the self-consistent mean field, 
the relativistic Hartree-Bogoliubov equations read 

liD-m-X A \ /Uk{r)\ /Uk{r)\ 

-A* -hj, + m + x) \Vk{r)) \Vk{r))' 

where hn is the single-nucleon Dirac Hamiltonian (^, and m is the nucleon mass. The 

chemical potential A has to be determined by the particle number subsidiary condition in 
order that the expectation value of the particle number operator in the ground state equals 



the number of nucleons. The column vectors denote the quasi-particle spinors and Ek are 
the quasi-particle energies. The pairing field A in (0) is defined 

Aa.(r, r') = W K.c.(r, r') ^ /7;,(r)\/,,(r'), (12) 

^ c,d Ek>0 

where a, b, c, d denote quantum numbers that specify the Dirac indices of the spinors, 
VabcdiT^-, r') are matrix elements of a general two-body pairing interaction. The RHB equa- 
tions are solved self-consistently, with potentials determined in the mean-field approxima- 
tion from solutions of Klein-Gordon equations for the meson fields. The current version of 



the model ||18[ describes axially symmetric deformed shapes. The Dirac-Hartree-Bogoliubov 
equations and the equations for the meson fields are solved by expanding the nucleon spinors 
t/fc(r) and Vfc(r), and the meson fields in terms of the eigenfunctions of a deformed axially 



symmetric oscillator potential [jT^. The calculations for the present analysis have been per- 



formed by an expansion in 12 oscillator shells for the fermion fields, and 20 shells for the 
boson fields. A simple blocking procedure is used in the calculation of odd-proton and/or 
odd-neutron systems. The blocking calculations are performed without breaking the time- 
reversal symmetry. 

The input parameters of the RHB model are the coupling constants and the masses for 
the effective mean-field Lagrangian, and the effective interaction in the pairing channel. In 
most applications we have used the NL3 effective interaction [|l9l for the RMF Lagrangian. 



Properties calculated with NL3 indicate that this is probably the best effective interaction 
so far, both for nuclei at and away from the line of /5-stability. For the pairing field we 
employ the pairing part of the Gogny interaction 



^ -({ri-r2)//i,) 
i=l,2 

with the set DIS for the parameters /Xj, Wi, Bi, Hi and Mj (i = 1,2). This force has 
been very carefully adjusted to the pairing properties of finite nuclei all over the periodic 
table. In particular, the basic advantage of the Gogny force is the finite range, which 
automatically guarantees a proper cut-off in momentum space. The RHB model with the 
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NL3+D1S effective interaction has been applied in studies of the neutron halo phenomenon 
in light nuclei []T5|, properties of light nuclei near the neutron-drip line ground state 



properties of Ni and Sn isotopes [E^, the deformation and shape coexistence phenomena 



that result from the suppression of the spherical N=28 shell gap in neutron-rich nuclei ^3 
the structure of proton-rich nuclei and the phenomenon of ground state proton emission 
18| , |2^p5| . In particular, it has been shown that neutron radii, calculated with the RHB 



NL3-I-D1S model, are in excellent agreement with experimental data |l21|j22 



III. NEUTRON DENSITIES FOR ATOMIC PNC EXPERIMENTS 

In this section we apply the relativistic Hartree-Bogoliubov model, with the NL3+D1S 
effective interaction, in the calculation of ground-state proton and neutron densities of Cs, 
Ba, Yb and Pb isotopes. These elements have been considered for measurements of isotope 
ratios of observables in atomic PNC experiments p|j5|,P-p^. 

Due to the large uncertainties in the experimental data on ground-state isovector prop- 
erties, most of the effective interactions used in the non-relativistic Hartre-Fock model, and 
the effective Lagrangians of the relativistic mean-field model, have not been specifically 
designed to describe neutron density distributions. Only more recently have isovector prop- 
erties been included in the set of data on which the effective interactions are adjusted. In 
the relativistic mean-field model, perhaps the most accurate set of meson masses and meson 
nucleon coupling constants is NL3 [T^. In Ref. this effective force has been used in 



the relativistic mean-field + BCS model to calculate the ground-state properties of 1315 
even-even nuclei with 10 < Z < 98. The calculated quantities include the total binding 
energies, rms radii, quadrupole and hexadecupole deformations. In order to illustrate the 
quality of the results that are obtained with the NL3 interaction, in Fig. |I] we display the 
binding energies and isotope shifts of the charge radii of Pb nuclei. The calculation has been 
performed with the RHB model, and the theoretical values are compared with the experi- 
mental data on binding energies and charge radii ^E\. The agreement between theory 
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and experiment is very good, and we notice in particular that the calculated radii reproduce 
the anomalous kink in the isotopic shifts at ^osp]^ [^,^]. The NL3 effective force has also 
been used in studies of properties of excited states. By using the time-dependent relativistic 
mean-field model and the relativistic random-phase approximation it has been 



shown that a self-consistent calculation with the NL3 interaction reproduces the excitation 
energies of giant resonances in spherical nuclei. Particularly relevant for our present study 
are isovector giant resonances. In Fig. ^ we show the isovector dipole and monopole strength 
distributions in ^^'^Pb, calculated with the NL3 effective interaction in the relativistic RPA 
framework. The calculated peak energy for the dipole resonance Ep = 12.9 MeV has to be 



compared with the experimental value of the excitation energy 13.3 ± 0.1 MeV ||3^. The 
calculated isovector monopole strength is much more fragmented, but the centroid is also in 
good agreement with the experimental value for the IV GMR: 26 ± 3 MeV |3^. In Fig ^ we 



plot the binding energies, quadrupole deformations, and isotope shifts of the charge radii for 
the Yb, Ba and Cs isotope chains. The results of the self-consistent RHB NL3+D1S calcu- 
lations are compared with the available experimental data P7|J2^J5^ . The overall agreement 



between the theoretical predictions and experimental data is excellent. The results shown 
in Figs. ^ demonstrate that the RHB model with the NL3+D1S effective interaction is 
adequate for the calculation of ground-state neutron densities. 

The self-consistent RHB neutron and proton ground-state density distributions for the 
isotopes 202,208, 2i4p|^ showu in Fig. ^. We will eventually compare the calculated rms 
radii for 202-2i4p]^^ notice the trend: the neutron radii increase with neutron 

number, while the neutron densities in the bulk region do not change much, except for some 
shell effects in the central region. The proton densities on the other hand, decrease in the 
bulk region as more neutrons are added on the surface. This effect, of course, has its origin 
in the symmetry energy. Therefore, one also observes an increase of proton radii, but much 
less pronounced than in the case of neutron densities. The difference in the trend of neutron 
and proton radii is emphasized in the inserts of Fig. ^ where the densities in the surface 
region are plotted on the logarithmic scale. 
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In Ref. [|T0| it was suggested that the wide range of stable isotopes of Yb might be used 
for a comparison of PNC observables. These nuclei are generally deformed and therefore, in 
order to represent graphically the deformed densities, the density functions are decomposed 



p{r) = ^Y.PL{r)YU9,<j,). (14) 



in multipole moments 



47r ^ 

The monopole L = 0, quadrupole L = 2, and hexadecupole L = 4 components of the neutron 
density in ^^^Yb are shown in Fig. ^ The deformations correspond to the neutron quadrupole 
moment Qn = 11.02 b, and the neutron hexadecupole moment if„ = 0.067 b"^. Since the 
hexadecupole deformations are generally small, in the following figures we compare only the 
monopole and quadrupole components of the ground-state proton and neutron densities. 
These are shown in Fig. ^ for the even-N Yb isotopes. The deformations generally decrease 
with increasing the neutron number (see the f32 values in Fig. The difference between the 
quadrupole components of the proton and neutron density is small for all Yb nuclei in the 
sequence = 98 — 106. The proton and neutron densities differ much more in the spherical 
part of the distribution, not only on the surface (different radii), but especially in the bulk 
region. Similar features are observed for the Ba (Fig. |^) and Cs (Fig. |]) isotopes. These 
nuclei are less deformed, and the spherical approximation for the densities will work better 
than in the case of the Yb isotopes. Correspondingly, the differences in the quadrupole 
components of the neutron and proton densities are smaller. The monopole components of 
these densities, however, display very different behavior both in the bulk and the surface 
region. 

The isotope shifts for the neutron radii of the ground-state density distributions of Cs, 
Ba, Yb and Pb nuclei are shown in Fig. |^. In all four cases a uniform increase of the 
neutron radii is observed. The differences between neutron and proton rms radii are plotted 



in Fig. |T0|, as function of the number of neutrons. The differences are of the order of 0.2 — 0.3 
fm and, of course, increase with the neutron number. The increase, however, is not very 
steep. For the Yb nuclei, for example, r^— increases by less than 0.05 fm between = 100 
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and = 106. The ratios r„/rp are compared in Fig. |TT]. We notice that while the isotope 
trend is similar, the actual values of r„/rp calculated with the RHB NL3+D1S model are 
considerably larger than those obtained in the Hartre-Fock Skyrme calculations for Pb |Q] 
and Cs 0,0. In particular, for ^^^Cs we calculate r^/r^ = 1.046, as compared to the Hartre- 
Fock values: 1.022 for SkM*, and 1.016 for Skill. An approximate estimate for the nuclear 
structure correction (^) to the weak charge has been derived in Ref. 

^Qw''^N{Zaf{.22lt)/q, (15) 

where the small parameter e is defined by 

{rl/rD ^ 1 + e. (16) 

The RHB NL3 value for r„/rp implies, therefore, that the nuclear structure correction to 
the weak charge can be more than a factor 2 larger than the estimate given in 0, bringing 
the value of AQi^T^ close to the uncertainties of atomic theory calculations. On the other 
hand, the results shown in Figs. 10 and ^ indicate that the changes in the neutron distri- 



bution parameters along the isotope chains should not present a limitation on the accuracy 
with which the Standard model can be tested in measurements of isotope ratios of PNC 
observables. 
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Figure Captions 

FIG. 1. Binding energies (left), and isotope shifts of the charge radii (right) of Pb nuclei, calcu- 
lated with the NL3 + Gogny DIS effective interaction. The theoretical values are compared with 



the experimental binding energies and charge radii |28| 



FIG. 2. The isovector dipole and monopole strength distributions in ^'^^Pb, calculated with the 
NL3 effective interaction in the relativistic RPA framework. 

FIG. 3. Binding energies (left), quadrupole deformations (center), and isotope shifts of the 
charge radii (right), for the Yb, Ba and Cs isotope chains. The results of the self-consistent RHB 



calculations are compared with the available experimental data 1 27, 34 25]. 



FIG. 4. Self-consistent RHB neutron and proton ground-state density distributions for 
202,208,2i4pj-j_ In the inserts the same densities in the surface region are plotted on the logarithmic 
scale. 

FIG. 5. The monopole, quadrupole and hexadecupole components of the self-consistent RHB 
neutron ground-state density distribution in the nucleus ^^^Yb. 

FIG. 6. Monopole and quadrupole components of the self-consistent RHB neutron (solid line) 
and proton (dot-dashed line) ground-state density distributions for the even-N Yb isotopes. 

FIG. 7. Same as in Fig. but for the Ba isotopes. 

FIG. 8. Same as in Fig. but for the Cs nuclei. 

FIG. 9. Isotope shifts for the neutron radii of the ground-state density distributions of Cs, Ba, 
Yb and Pb nuclei. 

FIG. 10. Differences between neutron and proton radii of the ground-state density distributions 
for the Cs, Ba, Yb and Pb isotope chains. 
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FIG. 11. Ratios of neutron to proton radii for the ground-state density distributions of the Cs, 
Ba, Yb and Pb isotope chains. 
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